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Abstract 

We have made an unbiased survey for molecular clouds in the Galactic Warp. 
This survey, covering an area of 56 square degrees at / = 252° to 266° and b = —5° to 
— 1°, has revealed 70 molecular clouds, while only 6 clouds were previously known in 
the region. The number of molecular clouds is, then, an order of magnitude greater 
than previously known in this sector at 14.5 kpc. The mass of the clouds is 
in a range from 7.8x10^Mq to 8.4x10^Mq, significantly less than the most massive 
giant molecular clouds in the inner disk, ~1O^M0, while the cloud mass spectrum 
characterized by a power law is basically similar to other parts of the Galaxy. The 
X factor, A^(H2)/ J T(^^CO)dV , derived from the molecular clouds in the Warp is 
estimated to be 3.5(±1.8) times larger than that in the inner disk. The total molecular 
mass in the Warp is estimated as 7.3 x 10^ Mq, and total mass in the far-outer Galaxy 
{R ^ 14.5 kpc) can be estimated as 2x10^ Mq. The spatial correlation between the 
CO and HI distribution appears fairly good, and the mass of the molecular gas is 
about 1% of that of the atomic gas in the far-outer Galaxy. This ratio is similar to 
that in the interarm but is ten times smaller than those of the spiral arms. Only 6 
of the 70 Warp clouds show signs of star formation at the IRAS sensitivity and star 
formation efficiency for high-mass stars in the Warp is found to be smaller than those 
in other molecular clouds in the Galaxy. 
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1. Introduction 



Molecular clouds are the sites of star formation, and play a crucial role in controlling 
the evolution of galaxies. It is therefore important to understand the physical conditions of 
molecular clouds and star formation therein under various physical conditions. 

Previous efforts to study molecular clouds have been mostly concentrated toward the 
inner Galaxy. These studies revealed the existence of giant molecular clouds (GMCs) having 
masses of 10^ to 10^ Mq, and that GMCs are actively forming both high- mass and low-mass 
stars (e.g., Mooney & Solomon 1988; Blitz 1993). The outer part of the Galactic disk is probably 
under significantly different physical conditions compared to the inner disk in several respects, 
including weaker stellar gravity, lower HI density, less intense UV fields and smaller cosmic 
ray flux (e.g., Bloemen et al. 1984). It is also known that the metallicity becomes lower in 
the outer Galaxy (e.g., Rudolph et al. 1997). Therefore, naturally we need to study the outer 
Galaxy, a unique nearby laboratory, to make new tests on how properties of star formation such 
as mass function, star formation efficiency etc. may change, depending on the above physical 
parameters. 

Optical components in the outer Galaxy are being revealed gradually. A deep photo- 
metric survey of late-type stars (mostly G and K dwarfs) toward the Galactic anticenter by 
Robin, Creze, & Mohan (1992) indicates a sharp cutoff of the old stellar disk beyond R = 14 
kpc. On the other hand. Population I stars apparently have a greater radial extent although 
they are not so many. Some HII regions have been found and studied at the outer edge of 
the Galaxy at R = 15-19 kpc, suggesting that star formation is taking place there at a lower 
level (e.g., Vilchez & Esteban 1996; de Geus et al. 1993; Kobayashi & Tokunaga 2000; Santos 
et al. 2000). It is thus very intriguing to study the dense molecular clouds and star formation 
in the outer Galaxy. 

To reveal the molecular distribution in the outer Galaxy, several efforts have been made 
to date to search for molecular clouds (e.g., Kutner & Mead 1981; Wouterloot & Brand 1989; 
Sodroski 1991; Brand & Wouterloot 1994; Digel, de Geus, & Thaddeus 1994; May, Alvarez, 
& Bronfman 1997). Kutner & Mead (1981) were the first to report the detection of low-level 
CO emission from molecular clouds in the first quadrant of the outer Galaxy. Although these 
studies cover a range of galactocentric distances (i?) from 8.5 kpc to 20 kpc and distances from 
the Galactic plane {\z\) of up to 1.5 kpc, they are limited in spatial resolution (May, Alvarez, 
& Bronfman 1997; Sodroski 1991) or spatial coverage (Digel, de Geus, & Thaddeus 1994) or 
are biased toward bright IRAS point sources (Wouterloot & Brand 1989; Brand & Wouterloot 
1994). Wouterloot & Brand (1989) observed ^^CO{J = 1-0) emission with SEST and IRAM 
telescopes in the direction of 1302 cold IRAS sources in the second and third Galactic quadrants 
and gave a catalog of those clouds (WB clouds, hereafter), providing a most comprehensive 
view of star-forming molecular clouds in the outer Galaxy (Wouterloot & Brand 1989; Brand 
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& Wouterloot 1994). We shall hereafter refer to the outer edge of the part, R ^ 14.5 kpc, as 
the far-outer Galaxy following Brand & Wouterloot (1995) who first used this word for R ^ 
16 kpc. Nonetheless, it is still not well known if there is a significant amount of non-star- 
forming molecular gas in the outer Galaxy. The importance of an unbiased study of molecular 
clouds covering non-star-forming clouds has been demonstrated by a number of authors over 
a large mass range, including those papers on cloud cores (e.g., Onishi et al. 1996; Onishi et 
al. 1998; Onishi et al. 2002) and on the giant molecular clouds in the Large Magellanic Cloud 
and M33 (Fukui et al. 1999; Engargiola et al. 2003). 

Most recently, Snell, Carpenter, & Heyer (2002) studied the properties of star-forming 
regions associated with 63 outer Galaxy [R > 13.5 kpc) molecular clouds, among which the 
most massive ones have molecular masses of ~10^ Mq, within a 60 deg^ area toward the second 
Galactic quadrant based on an unbiased CO Survey of the outer Galaxy at FCRAO (Heyer 
et al. 1998; Heyer, Carpenter, & Snell 2001). The survey has shown that the number of the 
molecular clouds dramatically decreases beyond R = 13.5 kpc. About 10 clouds among the 
63 have been found to be forming stellar clusters as revealed by K'-band imaging, and 39/246 
clouds at i? > 13.5 kpc have 2MASS extended source(s) that are candidates for star-forming 
sites (Snell, Carpenter, & Heyer 2002). Brunt, Kerton & Pomerleau (2003) decomposed the 
FCRAO data by applying a lower threshold for radiation temperature and doubled the sample 
of molecular clouds in the far-outer Galaxy in the second Galactic quadrant. These 2 papers 
show a possibility that a significant fraction of molecular clouds quiet in star formation still 
remains unobserved in the far-outer Galaxy. It is therefore an important task to better and 
firmly establish the properties of both star-forming and non-star-forming molecular clouds in 
order to have a comprehensive understanding of the molecular gas in the outer Galaxy. This 
goal should be attained by a large unbiased-sample of molecular clouds in the far-outer Galaxy. 

In this paper we describe the first results of the CO J = 1-0 survey of the third Galactic 
quadrant covering 56 square degrees for a velocity range corresponding to the far-outer Galaxy 
at i? = 14.5-19 kpc. We shall here focus on the southern Warp which shows a large shift from 
6 = 0° by keeping in mind that the Warp is potentially an ideal place to avoid contamination 
from the Galactic disk at every wavelength. The present work should be complementary to the 
work by Snell, Carpenter, & Heyer (2002) which studied the northern Warp. In this paper, 
details of the observations are given in section 2 and the results in section 3. Analysis and 
discussion are given in sections 4 and 5, respectively. Conclusions are given in section 6. 

2. Observations 

Observations were made in the J = 1-0 ^^CO emission line at 2.6 mm wavelength. The 
observed region is a part of the Warp in the outer edge of the Galaxy. The distribution of the 
Galactic Warp was first recognized in the 21 cm HI emission (Kerr 1957) and Burton & te Lintel 
Hekkert (1986) revealed a pair of significantly warped regions in the HI disk up to 1.5 kpc in 
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\z\. More recent papers on the HI Warp may be found elsewhere (e.g., Garcia-Ruiz, Sancisi, 
& Kuijken 2002). The present region is the southern part of the Warp, located mainly in the 
Galactic third quadrant. The Parkes 21 cm HI survey (Kerr et al. 1986) was used to select 
the region at / = 252°-266°, and h = —5° — 1°, where the HI emission is relatively strong at R 
> 14.5 kpc and the angular deviation of the Warp from 6 = 0° is largest. The area observed 
corresponds to about 56 square degrees in the sky. 

Observations were made with the NANTEN telescope at the Las Campanas Observatory 
in Chile in the period from July 2001 to August 2002. The main-dish diameter of the telescope 
is 4 m, providing a half-power beamwidth (HPBW) of 2.6 arcmin at 2.6 mm wavelength. This 
corresponds to a spatial resolution of 7.5 pc at a distance of 10 kpc, a typical distance to the 
southern Warp region. The front-end was an SIS mixer receiver. System temperature including 
the atmosphere was typically 200-250 K in SSB toward the zenith (Ogawa et al. 1990). We used 
two types of spectrometers. One of them was an acousto-optical spectrometer (AOS) which 
provided a velocity resolution of 0.65 km s~^ with a velocity coverage of 650 km s~^ (the Wide 
Band (WB) AOS). The other had a velocity resolution of 0.1 km s~^ for a velocity span of 100 
km s~^ (the Narrow Band (NB) AOS). The pointing accuracy was better than 20 arcsec, as 
checked by optical observations of stars with a CCD camera attached to the telescope, as well 
as by radio observations of Jupiter, Venus, and the edge of the Sun. In reducing the spectral 
data, we subtracted second-order or third-order polynomials determined from the emission-free 
velocity parts in order to ensure a flat spectral baseline. 

An absolute intensity calibration to obtain was made in the chopper-wheel method. 
In order to convert into we assume that target sources are extended compared with beam 
size. We observed Orion KL [a(1950) = 5^32"^47.^0, 5(1950) = -5°24'21"] every two hours, and 
the peak was taken as 65 K. All the observations were performed in the position-switching 
mode. "Off" positions, chosen to be as close as possible to "On" positions within two degrees, 
were checked to be emission-free at rms noise level of 0.07 K/ch in for both (WB and NB) 
AOSs. First, we observed the whole region at a 4 arcmin grid spacing with WB AOS. The 
coverage is thus under-sampled with a sampling fraction of 33%. The rms noise level in these 
observations was 0.2 K for 0.65 km s~^ resolution. Next, we observed the regions of CO emission 
detected at a 2a noise level at Vlsr > 90 km s~^ at a 2 arcmin grid spacing with NB AOS. 
The rms noise fluctuations were 0.25 K for 0.1 km s~^ resolution. If we adopt a fiat rotation 
curve (i?0 = 8.5 kpc, Vq = 220 km s~^, Kerr & Lynden-Bell 1986), the velocity threshold of 
90 km s~^ corresponds to R = 14.9 kpc at / = 252° and R = 14.4 kpc at / = 266°, respectively, 
and the maximum velocity covered Vlsr = 160 km s~^ corresponds to i? = 36 kpc at / = 252° 
and -R = 31 kpc at / = 266°, respectively. The total numbers of the observed positions in 
the 4 arcmin and 2 arcmin observations were ~13000 and ~7500, respectively. We admit that 
the under-sampling in the first step observations places a lower limit for the detectable cloud 
mass. Nonetheless, as is shown later in section 3.1, we are confident that the present survey is 
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complete in detecting molecular clouds whose mass is greater than 3000 Mq. 
3. Results 

3.1. CO distribution and cloud catalog 

Figure 1 shows the velocity integrated intensity distribution of CO in the Warp at Vi^^K 
> 90 km s~^. We have detected significant CO emission towards 702 positions at the 3a level, 
~0.9 K km s~^ in CO integrated intensity, and have identified 70 individual molecular clouds 
at the same level (0.9 K km s~^). The CO emitting area with the present beam is 0.78 deg^, 
1.4% of the observed area, indicating that the molecular clouds are very sparsely distributed 
in the region. Among the 70 clouds, 64 are newly detected, corresponding to 91% in number, 
and six of the 70 molecular clouds had already been detected and/or mapped by the previous 
observations. The total molecular mass of the 70 clouds is estimated to be 7.3x10^ Mq (see 
section 5.2), as compared to ~lxlO^ Mq previously known. In this work we have discovered 
molecular clouds an order of magnitude greater than previously known in this sector both in 
number and in mass, as explained in detail later in section 3.2. In Figure 2a we present some 
typical CO profiles which exhibit peak of 1-2 K with a small linewidth of a few km s~^. The 
present completeness limit in Leo is estimated as ~300 K km s~^ pc^, corresponding to ~3000 
Mq, which is set by the 4 arcmin grid and rms noise level of the first-step. On the other hand, 
in the second-step the present detection limit corresponds to Leo ~30 K km s~^ pc^ at L* = 
10 kpc and Leo ~70 K km s~^ pc^ at D = 15 kpc. The small cloud shown in Figure 2b is an 
example of a low mass cloud close to the detection limit in mass, ~800 Mq. 

Table 1 lists the observed properties in CO and physical parameters of the 70 Warp 
clouds as follows; (1) cloud number, (2-3) peak position in Galactic coordinates, (4) central 
velocity with respect to the local standard of rest (a composite profile is fitted by a Gaussian 
curve to derive this quantity), (5-6) distance from the Sun as determined from the rotation 
curve and from the Galactic center projected at 6 = 0° plane, (7) distance from the Galactic 
mid plane {b = 0°), (8) total integrated intensity, (9) full linewidth at half maximum (from 
Gaussian fit to a composite spectrum), (10) radius r as derived from r = (S'/vr)^/^ where 5* 
stands for the area at the three sigma intensity level, (11) beam-size corrected radius (here 
we assume that the beam pattern and cloud distributions are Gaussian, and cloud radii were 
corrected by using the following equation; (corrected size: Tc)^ = (observed size: r)^ — (beam 
size)^). (12) CO luminosity, Leo, calculated as (the total integrated intensity) x (grid area), 
(13) mass, Mco, derived from Leo by using the X factor (see section 4 for details on the X 
factor), (14) virial mass estimated as Mvir = 210x rc(size [pc])x dV^^ (linewidth [km s~^]) by 
assuming that the cloud is a uniform sphere of radius (beam-size corrected radius: Tc), and (15) 
comments; "cold IRAS" indicates clouds which are associated with a cold IRAS point source 
(selection criteria; F25 < Fqq with qualities > 2 at wavelengths of 25 and 60 micron) near the 
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peak. 

3.2. Comparison with the previous works 

We shall compare the present CO results with the previous works. Brand & 
Wouterloot (1994), BW94 hereafter, selected and mapped twenty-seven WB clouds at 
16 kpc, including thirty star-forming IRAS point sources. Among those, Wouterloot & 
Brand (1989) detected 6 molecular clouds at Vlsr > 90 km s~^ in the present observed re- 
gion (WB1096, 1103, 1121, 1126, 1152, 1153) and BW94 mapped four of them in CO J = 
1-0 emission. The present work has detected four of the 6 clouds and these coincide with the 
WB clouds, respectively, as follows; No. 21 (WB1096), No. 27 (WB1126), No. 63 (WB1152) 
and No. 64 (WB1153). No. 71 (WB1077), located on the border of the present area, has not 
been fully mapped. The present table does not include the remaining two clouds (WB1103 
and WB1121). WB1103 was mapped by BW94 and its Leo is calculated to be 5.1 K km s~^ 
pc^, below the present detection limit. WB1121 was observed only toward one point, the IRAS 
source position, by Wouterloot & Brand (1989) and was not mapped by BW94. The single 
point CO spectrum gives a lower limit of Leo as ~30 K km s~^ pc^, somewhat smaller than 
the present detection limit, and WB1121 is perhaps missed because of its low Leo- 

Our observations do not spatially resolve small structures of less than ~10pc. Toward 
four molecular clouds No. 21, 27, 63 and 64, Wouterloot & Brand (1989) made CO observations 
with the 15 m SEST telescope and gave peak T^, (roughly equal to T^,) equal to 2.2-4.5 K, 
which is a few times higher than we obtained. This indicates that the small-scale distribution 
within the clouds is beam-diluted. On the other hand, a lower resolution study at 8.8 arc-min 
beam by May, Alvarez, & Bronfman (1997) detected only one cloud in this area at / = 257.5°, 
having 1.5x10^ Mq, corresponding a present cloud complex resolved into three clouds. No. 27, 
28 and 30, with the 2.6 arcmin beam. 

3.3. Mass spectrum 

The cloud mass, Mqo, ranges from 7.8x10^ Mq to 8.4x10"^ Mq, and the largest mass 
is significantly smaller than that of giant molecular clouds in the inner disk, i.e., ~10^ Mq. 

The mass distribution of 68 clouds out of 70 is shown in Figure 3. Two clouds (No. 
18 and 56) are not included here because their masses are not well determined due to double 
peaked profiles (see also Table 1). The mass spectrum of the 47 clouds is well fitted by a 
single power law for a mass greater than 3000 Mq by taking the maximum-likelihood method 
(Crawford & Jauncey 1970), N(> Mco/Mq) = 1.3x10^ Mc5-5.8, where a = -0.68±0.16. 
The mass spectra studied for the Galactic CO clouds are well fitted by a power law with an 
index value of a ~0.5-0.8 (e.g., Heyer, Carpenter, & Snell 2001; Yonekura et al. 1997; Digel et 
al. 1996; Brand & Wouterloot 1995; Solomon & Rivolo 1989). The present result is consistent 
with these studies within the error. This leads to the conclusion that the cloud mass spectrum 
characterized by a power law is basically the same throughout the Galaxy. The slope of the 
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mass spectrum for clouds with < 3000 Mq becomes shallower, indicating that the sampling is 
incomplete below this mass. 

S.Jf.. Overall uncertainties of the cloud mass 

Individual Warp clouds show Vlsr around 100 km s~^, and their galactocentric distances 
are estimated to be i? = 14.4-18.6 kpc. The distance given in Table 1 may bear errors by the 
following two reasons, the uncertainties in the rotation curve and the cloud-to-cloud random 
motion. An assumption of a flat rotation curve may result in error up to about 15% in distance 
(see e.g.. Brand & Blitz 1993). We note here that the rotation curve in the far-outer Galaxy is 
only poorly determined at present and that a fiat rotation curve is fairly commonly observed 
in the outskirts of external spiral galaxies (see for a review e.g., Sofue & Rubin 2001). We 
also note that the present targeted region is less affected by uncertainties in kinematic distance 
compared to the regions close to the Galactic center or to the anti-center where the rotation 
curve affects the distance estimate more significantly. The other source of error is the random 
motion of the clouds. If we adopt 5 km s~^ as the typical cloud-to-cloud velocity dispersion as in 
the inner disk (e.g., Israel et al. 1984), the distance ambiguity may be up to ~8% of those given 
in Table 1. Another possibility for uncertainties in distance is systematic large non-circular 
motions such as a streaming. This is difficult to estimate quantitatively, and is assumed here to 
be negligible because streaming is usually prominent only in the vicinity of spiral arms. Thus, 
the overall uncertainty is up to 17% in distance. Another source of uncertainty is calibration 
of the absolute fiux of CO, about 10%. By including all these, the total uncertainty in mass is 
estimated to be < 38% for a given X factor. 

4. Derivation of the X factor 

We shall estimate the X factor (the conversion factor), X = N{'H2) / 1 T{'^^CO)dV 
[cm~^/K km s~^], in the Warp. This is practically a useful factor to convert the CO luminosity 
to molecular mass and is also a fundamental factor that may reflect the physical properties of 
molecular clouds (e.g., Paglione et al. 2001; Strong et al. 1988). The X factor for the Galactic 
disk clouds is around 2x10^° [cm~^/K km s~^]. The X factor is also written as X = Mqq/ Lqq 
[Mq/K km s-i pc2], e.g., Mco = 4.3xLco (if AT = 2x102° [cm'^/K km s"!]). 

We show an Lco-M^h diagram in Figure 4. We have plotted 54 Warp clouds that have 
more than three consecutive observed points in order to reduce the effects of beam dilution in 
cloud size. In addition, we have included, for comparison, 216 Carina arm clouds (see section 
5.1) observed with the same instrument (Matsunaga et al. 2005, in preparation) and 75 GMCs 
in the inner disk with known distances (see Table 1 in Solomon et al. 1987 and Appendix). In 
this figure, we used size-corrected clouds (see section 3.1) both for the Warp and the Carina 
regions. All virial masses are derived from the calculation Mvir = 210xrxdV'^ by assuming 
for simplicity that density distribution is uniform in a cloud. The Warp clouds are indicated 
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by crosses, those in the Carina arm by open circles and those in the inner disk by diamonds. 
For all the three regions, Leo and Mvir show a good positive correlation with a slope of ~1. 
This figure also indicates that Leo is significantly different among the three regions for a given 
Myir- The Warp clouds generally show Leo a factor of two smaller than that of the Carina 
clouds and a factor of three smaller than that of the inner disk clouds; these differences may 
depend on the range in Mvir- It is to be noted that the slopes of the Lco-^vir relation for the 
three regions are slightly different, 1.05(±0.13) for the Warp region, 0.93(±0.04) for the Carina 
region and 0.78(±0.03) for the inner disk. 

We shall estimate the difference in the X factor by comparing the offset in Leo foi" a 
given virial mass. We restrict ourselves to the mass range 10^'°-10^'^ Mq by keeping in mind 
that the accuracy of mass for smaller clouds whose Mvir is less than 10^ Mq is worse than 
clouds whose M^ir is greater than 10^ Mq. In this range, there are still a significant number 
of sample clouds (see Figure 4), i.e., 33 Warp clouds and 20 inner disk clouds, allowing us to 
make a reasonable comparison with statistical significance. The solid lines in Figure 4 indicate 
regression lines with a slope of 1.0 for the Warp and the inner disk clouds; they are expressed 
as log(Mvir) = log(Leo) + l-41(±0.25) and log(Mvir) = log(Leo)+0.86(±0.24), for the Warp 
and the inner disk, respectively. The averaged value of Leo is then estimated to be a factor 
3.5±1.8(1 0") smaller than that for the inner disk. 

In this analysis we did not make a correction for size for the inner clouds, nor take into 
account the possible overestimation of linewidth due to cloud overlap in the line-of-sight for all 
regions. If we consider these effects, the virial masses of inner clouds would decrease, resulting 
in a larger increase in the X factor of the Warp. 

Gamma-ray observations have been used to estimate the X factor in the inner disk, X 
= 1.6x1020 [cm-VK km s'^] (Hunter et al. 1997) and X = 1.9x10^° [cm'^/K km s'^] (Strong 
& Mattox 1996). We shall adopt here the former value of the X for the inner disk, 1.6x10^° 
[cm~2/K km s~^], and the above estimate of the difference in Leo, 3.5. By taking these values 
we obtain the present X = 5.6 xlO^'' [cm~2/K km s~^] in the Warp. 

It is suggested that the X factor in the outer Galaxy (9 kpc ^ i? ^ 14 kpc) is a few times 
larger than in the inner Galaxy (Mead & Kutner 1988; Sodroski 1991; Digel et al. 1996). The 
present X factor at 14.5 kpc is consistent with these, indicating that X factor becomes 
larger by a factor of 2-4 at i? ^ 9 kpc. The difference of the X factor may be understood as 
due to lower temperature (see next Sect. 5.1) and/or to lower metallicity in the Warp. In this 
paper, we calculate cloud masses by using X = 5.6 xlO^^ [cm~^/K km s~^] or the relationship 
of Meo [Mq] = 12.1xLeo [K km s"^ pc^]. 

The present derivation of the X factor is based on the assumption that virial mass (M^ii.) 
remains the same for a cloud of given mass. Almost all the CO clouds in the Warp have larger 
Mvir compared with Meo and the ratio of Myir/Meo ranges from 0.5 to 6.4 with an average 
value of 2.1 (Table 1). It is possible that the clouds are in dynamical equilibrium under the 
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external pressure or that they are transient features with a timescale of the order of the typical 
crossing time, 3x10^ yrs. Smaller clouds with Leo < 3x10^ K km s~^ pc^, in particular, may 
not be bound by self-gravity (Heyer, Carpenter, & Snell 2001). The virial mass range used in 
the present work > 10^ Mq perhaps guarantees the above assumption. 

The present result on the X factor is consistent with previous work by Mead & 
Kutner (1988); Sodroski (1991); Digel et al. (1996) for the outer Galaxy. On the other hand. 
Brand & Wouterloot (1995), BW95, and Snell, Carpenter, & Heyer (2002) derived a conclusion 
that the X factor is not significantly different between the far-outer Galaxy and the inner disk. 
These authors used ^^CO emission to estimate the molecular hydrogen column density by using 
an LTE analysis with an H2/^'^C0 ratio similar to that in the inner disk (Snell, Carpenter, 
& Heyer 2002) and with a higher H2/^^C0 ratio extrapolated from ^^CO integrated intensity 
(Brand & Wouterloot 1995). They argued that the thus derived molecular column density 
agrees with that derived from ^^CO integrated intensity using the conventional X factor of 
the inner disk, ~2xl0^° [cm~^/K km s~^], and concluded that the X factor is uniform in the 
Galaxy. This argument however may not be justified because the difference in the density 
regime probed in ^^CO and ^^CO is not properly taken into account. It is well known that 
the ^^CO emitting region is larger than the ^^CO emitting region in a given cloud; a typical 
mass ratio of ^^CO and ^^CO emitting regions are ~3 in the local clouds (see e.g., Mizuno et 
al. 1998; Mizuno et al. 2001; Onishi et al. 1999; Tachihara et al. 2001; Hara et al. 1999), and 
this is generally interpreted in terms of different self-shielding against far ultraviolet radiation 
(e.g., Warin, Benayoun, & Viala 1996). In the velocity space, this difference manifests itself as 
the broader line profile in ^^CO than in ^^CO, and the ^^CO line profiles indeed exhibit emission 
in a velocity range where no ^^CO emission is seen. Taking this into consideration, we need to 
accept that the molecular column density derived from ^^CO refers to only the ^^CO emitting 
region, although the whole ^^CO emitting region should be more extended and more massive 
than the ^^CO emitting region. Therefore, the molecular column density derived from ^^CO 
should be smaller than that derived from ^^CO reflecting the ^^CO/^'^CO cloud mass ratio, ~3; 
it is likely that the exact ratio in column density may depend on the cloud geometry. In order 
to reconcile this with the results of BW95 and Snell, Carpenter, & Heyer (2002), we have two 
alternatives; one is that the ratio Ha/^^CO may be smaller than BW95 assumed and the other 
is that the X factor is larger in the Warp than in the inner disk. The latter is actually in 
accord with the present result and is perhaps more likely (see also Sect. 5.1). We thus confirm 
the present conclusion that the X factor in the Warp is a factor of 3.5 larger than in the inner 
disk, noting that the conclusion by BW95 and Snell, Carpenter, & Heyer (2002) needs to be 
carefully reexamined by considering the ^^CO/^^CO cloud mass ratio. 
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5. Discussion 



5.1. Comparison of the Warp clouds with the Carina- Arm clouds 

In order to compare the Warp clouds with the clouds in the Galactic disk, we have 
created a diagram of physical parameters that include clouds in the Warp and the Carina arm 
as shown in Figures 5. We have chosen here the Carina arm because the region was observed 
with the same instrument, NANTEN (Matsunaga et al. 2005, in preparation), and because 
cloud parameters can be better defined with less confusion than in the inner disk. The CO data 
of Carina were taken at rms sensitivity ~0.3 K with WB AOS, 3 times worse than that for the 
Warp clouds, and at a 4 arcmin grid spacing, twice of that for the Warp. These observational 
parameters place the lower limit in cloud mass 20 times higher than in the Warp. We have 
thus used only the reliable sample clouds in Carina whose mass is greater than 16,000 Mq and 
restricted the cloud distances from the Sun {D) in a range from 10.3 kpc to 15.3 kpc, the same 
range for the Warp clouds. 

Figure 5 gives the size-linewidth relation for the same sample of CO, a relationship often 
used to characterize the physical properties of clouds. In this figure, we also plotted the beam- 
corrected size (see Sect. 3.1) of the Warp clouds (crosses) and the Carina arm clouds (circles) 
in addition to the far-outer Galaxy clouds in the second Galactic quadrant observed by Heyer, 
Carpenter, & Snell (2001). The broken line and dotted line are least-squares fits calculated 
by BW95 for the inner- {R < 8.5 kpc) and outer-Galaxy clouds {R > 8.5 kpc), respectively. 
The diagram shows that molecular clouds in the far-outer Galaxy, from NANTEN and FCRAO 
data, are roughly consistent with the outer Galaxy line by BW95. This also clearly indicates 
that the outer Galaxy clouds are characterized by linewidths a factor of two smaller than the 
inner Galaxy clouds (see also Figure 10). On the other hand, the Carina clouds are located 
between the two lines for the outer Galaxy and for the inner Galaxy. If correct, this suggests 
that the clouds in the disk may have a gradient in linewidth as a function of R, which should 
be tested by a more extensive dataset in the future. 

Here, we will again discuss on the X factor based on the above relations in order to 
obtain a deeper insight into the cause of the larger X factor in the Warp. If we denote the 
cloud size in the far-outer Galaxy and in the inner disk as r^'^'^ and r^^, respectively, the 
size-linewidth relation in Figure 5 is expressed as d^ = A;/2x (r*^^^)^/^ and dV = kx(r^^Y^'^, 
respectively. On the other hand, BW95 (Fig. 13) showed that Leo oc dV^; we have confirmed 
this relation holds also for the present Warp clouds(see Appendix and Figure 10). By using the 
above equations and Myir oc rxdV'^, we find that Leo in the far-outer Galaxy is one fourth of 
Leo in the inner disk. This offers a simple explanation of how the X factor becomes larger by 
a factor of 3.5(~4) in the Warp, providing further support for the present X factor. 
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5.2. Comparison with HI 

The CO distribution is superposed on the HI (Kerr et al. 1986) distribution in the l-b 
diagram (Figure 6). This figure indicates that the HI and CO are well correlated with each 
other in the sense that the CO clouds, while much clumpier than HI, are distributed mostly 
toward the intense HI emission. To be more quantitative, 93% (= 65/70) of the CO clouds are 
located within the 800 K km s"^ contour level of HI in Figure 6. 

Figure 7 shows the distance z from the Galactic mid plane as a function of R. The peak 
positions of CO clouds are superimposed on an HI intensity map integrated between / = 252° 
and 266° from the HI data. This clearly shows that the CO clouds are distributed up to i? = 
18.6 kpc and there is no molecular cloud beyond this radius as is consistent with BW94. The 
trend that \z\ increases with R is clearly seen, an obvious sign of the Warp. Wouterloot et 
al. (1990) also showed the sign of warp structure by CO and HI observations on an l-z map. 
Comparisons between the HI and CO in Figures 6 and 7 indicate that the Warp clouds are 
not distributed uniformly, but that they are distributed in some groups extended over a few 
degrees, coinciding with larger HI distributions of similar extents. In Figure 7 one also sees 
that the ^-distribution of the molecular clouds becomes wider with increasing R, a sign of the 
flaring of the molecular disk. 

Figure 8 shows a face-on view of the outer disk. The peak positions of CO clouds in the 
Warp and the Carina arm are superimposed on an H I column density map from the H I data at 
^LSR > 10 km s~^ (Kerr et al. 1986). We used an intensity-mass conversion factor of 1.82x10^^ 
[cm"^/K km s~^] to calculate HI mass. Thin broken contours indicate a column density of 
0.5x10^^ cm~^. For the Warp and the Carina arm, the total molecular mass is estimated as 
7.3x10^ Mq and 2.1x10^^ M©, respectively, while that of HI is calculated to be 4.4x10^^ M© 
and 14.5x10^ M©, respectively. We applied a standard conversion factor, X = 2x10^° [cm~^/K 
km s~^], for the Carina clouds to calculate their masses. By using the total area (11 kpc^ for 
the Warp and 24 kpc^ for the Carina arm), averaged surface densities of H2 and HI in the Warp 
and the Carina arm are calculated within the thick broken line and are listed in Table 2 along 
with the corresponding values in other regions of the Galaxy, i.e., the inner disk, the Perseus 
arm and the interarm, and in the LMC for reference. 

We find that the molecular cloud mass in the Warp is ~1% of the total hydrogen mass. 
This ratio is significantly less than that in the Carina arm (14%) or in the Perseus arm (11%) 
or in the inner disk (50%), but is similar to that in the interarm (1%). The ratio is also much 
less than that in one of the nearest galaxies, the LMC (Mizuno et al. 2001). This result is 
consistent with the conclusion reached by Snell, Carpenter, & Heyer (2002) for the northern 
far-outer Galaxy and suggests that the conversion efficiency from atomic to molecular gas is ten 
times smaller in the far-outer Galaxy than in the inner disk. The disk layer thickness (FWHM) 
of HI gas increases with R, 400 pc at i? = 10 kpc and 600 pc at i? = 16 kpc (Wouterloot et 
al. 1990), and the HI mass density is calculated as 1.5x10'' M©/kpc^ and 0.7x10''' M©/kpc^ 
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for the Carina arm and the Warp, respectively. We suggest that this HI density decrease by a 
factor of 2 may lead to the lower efficiency of conversion into molecular gas in the Warp. 

We shall discuss the total molecular mass in the far-outer Galaxy. The total molecular 
mass is estimated to be 7.3x10^ Mq in the present region of the southern Warp. The area 
of the present survey corresponds to 3.5% of the area at i? = 14.5-19 kpc. If we assume the 
proportion of molecular cloud in the far-outer Galaxy is uniform, the total molecular mass in 
the far-outer Galaxy is estimated to be 2.1x10^ Mq by simply multiplying 100/3.5 to 7.3x10^ 
Mq. We can also check the total molecular mass in the same way by using the FCRAO data 
from the on-line catalog (Heyer, Carpenter, & Snell 2001; Brunt, Kerton & Pomerleau 2003). 
Heyer, Carpenter, & Snell (2001) found 117 clouds at i? > 14.5 kpc within a region extending 
40 deg in longitude, three times larger than the area in our survey, and gave Leo for each. The 
area covered corresponds to ~6% of the far-outer Galaxy. The total molecular mass in that 
area is calculated to be 5.0x10^ Mq (total CO luminosity is 4.1x10^ K km s~^ pc^), if the 
same X factor as in the Warp is applied. By assuming a uniform molecular cloud distribution 
over the far-outer Galaxy, the total molecular mass is estimated to be 8x10® Mq, somewhat 
smaller than the above for the southern Warp. This may be due to a moderate asymmetry 
in molecular mass between the south and north, while the difference is only within a factor of 
two. We thus infer that the current best estimate for the total molecular mass in the far-outer 
Galaxy is 1-2x10^ Mq. The total HI mass in the far-outer Galaxy is ~2xl0^ Mq as derived 
by a few groups (Wouterloot et al. 1990; Wolfire et al. 2003), and the overall molecular /atomic 
ratio is 0.5-1%. On the other hand, Bronfman et al. (2000) estimated the total molecular mass 
at 8.5 <R< 14.5 to be 2x10^ Mq (X = 1.56x 10^° [cm^VK km s^^]). We thus infer that 
the total molecular mass in the far-outer Galaxy is 5-10% (depending on X; see Bronfman et 
al. 2000) of that in the outer Galaxy. 

5. 3. Star formation in the Warp clouds 

In order to investigate star formation in the Warp clouds we have used the IRAS point 
source catalog to select possible candidates for YSOs associated with the 70 Warp clouds. 
Cold IRAS point sources (F25 < Fqq with qualities > 2 at both bands 25 fim and 60 /im) are 
associated with the peak CO emission in only 6 clouds (see Figure 1). These clouds may be 
forming massive to intermediate mass (B type) stars as judged from their IRAS flux density; we 
roughly estimate the total far-infrared luminosity of these IRAS sources from the flux densities 
at 12-100 fim, as Lpjji = 2-19x10^ Lq, where we use a simple method of estimating Lpj^i from 
the sum of the fluxes at 4 bands (Casoli et al. 1986). 

The LpiFt/ Mqo ratio for these 6 clouds ranges from 0.1 Lq/ Mq to 2.7 Lq/ Mq with an 
average value of 0.8 Lq/ Mq. We note that this ratio is smaller than the previous studies. Snell, 
Carpenter, & Heyer (2002) estimated the L/M ratios, 0.4-33 Lq/Mq with an average value of 
8.1 Lq/Mq, based on 9 clouds in the far-outer Galaxy at i? = 13.5-17.4 kpc, where they used 
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an X factor of 1.9x10 , a factor of 2.9 smaller than the present value. If the same X factor 
for the Warp, X = 5.6x10^° [cm~^/K km s"^], is applied, averaged L/M decreases to be 2.8 
Lq/Mq. Snell, Carpenter, & Heyer (2002) suggested that the L/M ratios for the far-outer 
Galaxy clouds are similar to those for the inner disk clouds, which are in the range 0.5-30 
Lq/Mq with an average value of 4 Lq/Mq (Mooney & Solomon 1988). The present sample 
indicates that the L/M ratio in the southern Warp is a factor of 4-5 smaller than that in the 
northern far-outer Galaxy and inner Galaxy. We therefore suggest that the star formation 
efficiency in the southern Warp may tend to be smaller although we definitely need to pursue 
this issue by increasing the number of sample clouds with more sensitive infrared source list in 
the far-outer Galaxy. 

We have created a histogram of cloud mass. Figure 9 indicates that massive clouds 
having larger molecular column density form stars (dark gray) more often than lower mass 
ones, a similar trend to the disk cloud (e.g., Yonekura et al. 1997; Kawamura et al. 1998; Hara 
et al. 1999; Tachihara et al. 2002; Onishi et al. 2002). The trend seems to be consistent but not 
so strong as in the case of the disk cloud core, possibly due to different range of cloud mass, the 
lower spatial resolution of the present molecular data and the lower sensitivity of the infrared 
sources. We suggest that the present molecular clouds provide a set of good targets which can 
be used to make detailed searches for star formation at various wavelengths, in particular at 
near- to mid-infrared wavelengths. 

It has been shown by several authors that the ratio between the virial mass and the 
luminous mass of molecular clouds is a good indicator of the gravitational relaxation of clouds 
leading to star formation on large scale (e.g., Yonekura et al. 1997; Kawamura et al. 1998; 
Tachihara et al. 2002) and on small scale (e.g., Hara et al. 1999; Onishi et al. 2002). These 
studies show that the virial mass derived from the velocity dispersion of molecular spectra 
tends to become smaller for clouds with star formation, and that star formation takes place 
for clouds having greater molecular column density. These are reasonable characteristics that 
indicate relaxation of the cloud structure under self gravity. 

In order to test how the above ratio depends on star formation and cloud mass in the 
far-outer Galaxy, we have created a histogram of cloud mass along with the cloud mass ratios. 
Figure 9 clearly indicates that massive clouds having larger molecular column density form 
stars (dark gray) more often than lower mass ones, a similar trend to the disk clouds (see the 
references above). It is also seen that the virial to luminous mass ratio shown by squares and 
solid lines for the star-forming clouds tend to be somewhat smaller than that for non-star- 
forming ones (circles and broken lines). The trend seems to be consistent but not so strong 
as in the case of the disk clouds, possibly due to the lower spatial resolution of the present 
molecular data and the lower sensitivity of the infrared sources. We suggest that the present 
molecular clouds provide a set of good targets which can be used to make detailed searches for 
star formation at various wavelengths, in particular at near- to mid-infrared wavelengths. 
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6. Conclusions 



We summarize the present conclusions as follows; 

1) A first unbiased study of the molecular clouds in the southern Galactic Warp has 
detected 70 molecular clouds, nearly ten times more than previously known, in an area of 56 
square degrees, thereby providing more support to the conclusion that the molecular mass in 
the far-outer Galaxy (we defined as i? ^ 14.5 kpc) amounts to 2x10'' Mq. 

2) The conversion factor, A^(H2)/ J T{^'^CO)dV, in the far-outer Galaxy derived from 
Lco-^vir relation is 5.6x10^° [cm~^/K km s~^], 3.5(±1.8) times larger than that in the inner 
disk. 

3) The size-linewidth relation shows offset between the inner disk and far-outer Galaxy 
clouds in the sense that cloud size in the far-outer Galaxy is twice that in the inner disk at any 
given linewidth. 

4) The CO and HI show a good spatial and kinematical correlation although the CO 
clouds are much clumpier than H I. Warp structure and flaring are clearly seen in the molecular 
cloud distribution, as well as in H I. The total molecular mass in the far-outer Galaxy is about 
1% of that of HI, ten times lower than in the spiral arms and is similar to that in the interarm. 

5) Only 6 of the 70 Warp clouds show signs of star formation, i.e., far- infrared compact 
sources detected with IRAS, although the far- infrared observations are sensitive enough only 
for high mass stars of B type or of earlier. The star formation efficiency derived from this is a 
factor of 4-5 smaller than that in the northern counterpart of the Warp. 

Appendix. Comparison of CO data in the inner disk 

We have used the CO data of the inner disk clouds from the catalog of Solomon et 
al. (1987) to derive the X factor in the Warp (section 4). There is however a concern about the 
reliability of the Solomon et al.'s data because the clouds were defined at a certain high level of 
the temperature, i.e., 4 K (Solomon et al. 1987). This was unavoidable due to the continuous 
high intensity level of the CO emission in the inner disk, even in between the clouds. On the 
other hand, the present Warp clouds are distributed more sparsely, allowing a cloud definition 
at a very low CO integrated intensity of 0.9 K km s~^. A natural question here is if the rather 
high threshold value in the CO intensity in the Solomon et al.'s data may cause some errors in 
the cloud physical parameters in deriving the X factor. It is therefore important to test if it 
does not cause any systematic error to use the Solomon et al.'s results in deriving the X factor. 

We have used the NANTEN ^^CO Galactic Plan Survey results (/ < 60°) on the Solomon 
et al.'s clouds. From the NANTEN dataset we have defined 55 clouds which are coincident 
with the distance calibrator clouds in the Solomon et al.'s catalog which includes 273 clouds, 
and provided three diagrams in Figure 10, i.e., Lco-M^h, size-linewidth, and Lco-linewidth 
relations. Details of the definition of size and linewidth are not given in BW95, while the same 
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trend that the inner Galaxy clouds have large linewidths is clearly seen. For 11 clouds, we 
calculated Leo and Mvir at three different threshold values of CO integrated intensity in order 
to investigate the dependence on the observational sensitivity. Spectra of the NANTEN toward 
No. 181 cloud (from the Solomon et al's catalog), for instance, are shown in Figure 11, for the 
three lowest integrated intensity (and temperature) values of 53.8(4.6), 12.8(1.3), and 5.3(0.4) 
K km s~^(K), respectively. Figure 12 indicates the correlations between Leo and Mvir, showing 
that plots move along the least-squares fit for the 55 clouds, log(Mvir) = 0.91xlog(Lco)+l-ll 
(c.c = 0.91). 

This comparison demonstrates that the regression does not change depending on the 
threshold CO intensity among the NANTEN datasets and that the present analysis is fairly 
consistent with the Solomon et al.'s regression employed in the present work. The individual 
values of Leo and Mvir of course depend on the threshold and changes a little as shown in 
Figure 12. Keeping in mind that the determination of the X factor actually mainly depends 
on relative shift among the regressions in Figure 12, we conclude that the high threshold value 
in the Solomon et al.'s data are not likely the cause of any significant error in the present work. 
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Fig. 1. Velocity integrated intensity distribution of CO (J — 1-0) in the Warp region at I^lsr > 90 km 
s^^ (corresponding to i? = 14.9 kpc at I ~ 252° and R = 14.4 kpc at I = 266°, respectively). Contours 
are every 3.0 K km s^^ from 0.9 K km (Sc). Crosses indicate the position of cold IRAS sources that 
are associated with CO clouds. 
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Fig. 2. Distribution of CO integrated intensity and CO profiles of selected Warp clouds. Contours are 
every 1.2 K km s^^ from 0.9 K km s^^. The filled circles indicate the observed positions at 2' grid and 
cloud numbers from Table 1 are given in the individual panels. These figures show (a) typical clouds and 
(b) one of the smallest clouds in the present sample. 
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Fig. 3. Mass spectrum of 68 Warp clouds. The number of clouds, N(> Mco/Mq) with mass greater 
than Mco, is plotted against the cloud mass Mqo, along with the best-fitted power law (dashed line). 
The power law N(> Mco/Mq) = 1.3x10'' Mc^*^*=^°-^^-5.8 is derived by using the maximum likcUhood 
method of Crawford & Jauncey (1970). 
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Fig. 4. Virial mass as a function of CO luminosity. The Warp clouds are indicated by crosses. Gray 
open circles show the Carina arm clouds (see Figure 8), and diamonds show the inner disk 'calibra- 
tor clouds' with known distances (Solomon et al. 1987). Three interrupted lines are fitted by the 
least-squares method; Warp: log(A/vir) = 1.05(±0.13) xlog(ico) + l-19(±0.38), correlation coeflicicnt(c.c.) 
= 0.84; Carina: log(Mvir) = 0.93(±0.04)xlog(Lco)+l-28(±0.18), c.c. = 0.90; Inner-disk: log(Afvir) = 
0.78(±0.03)xlog(Lco)+l-70(±0.17), c.c. = 0.97. Two thick solid lines indicate regression lines with a 
slope of 1.0 for the Warp and the inner disk clouds in the mass range 10"' *'-10^'^ Mq; they are expressed 
as log(A/vir) = log(Lco)+l-41(±0.25) and log(Afvir) = log(Lco)+0.86(±0.24), respectively. 
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Fig. 5. Size-linewidth relation of the Warp clouds (crosses). For comparison, CO clouds in the Carina 
arm (far-side) at 9 kpc < i? < 12 kpc (circles) and from the FCRAO data with R > 14.5 kpc by Heyer, 
Carpenter, & Snell (2001) (diamonds) are shown. Broken and dotted lines in this figure for 'inner Galaxy' 
and 'outer Galaxy' are cited from Brand & Wouterloot (1995). 
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Fig. 6. The CO distribution is superposed on the H I distribution in Galactic coordinates. The gray-scale 
map indicates the HI velocity integrated intensity of 90-140 km s^^ in Vlsr (approximately corresponding 
to R ^^14.5-25 kpc), and thin dotted contour indicates the integrated intensity = 800 K km s^^ contour 
level of HI. Crosses indicate the positions of CO clouds. The broken line shows the observed area. 
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Fig. 7. The Warp clouds are shown in a plane; the distance from the Galactic center (R) vs. distance 
fr-om the plane of & = 0° (z). Crosses indicate the peak positions of CO clouds in the Warp (Vlsr > 90 
km s^^) and the gray-scale indicate HI distribution integrated over Galactic longitude from 252° to 266°. 
Contours show 30%, 45% and 60% levels of the maximum HI intensity (located at i? = 11.5 kpc and z — 
—250 pc). The observed area is denoted by solid lines. 
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Fig. 9. Mvir histogram for the Warp clouds with (dark gray) and without (hght gray) IRAS sources (left 
axis). Averaged Mvir/Mco ratio for clouds with (squares and solid lines) and without IRAS (circles and 
dotted lines) in each bin are presented at the right axis. 
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Fig. 11. Spectra from the NANTEN data toward No. 181 cloud, which number is cited from Table 1 of 
Solomon et al. (1987). X mark represents the center position which Solomon et al. (1987) gave. Dotted 
lines show the center velocity, Vlsr = 97 km s^^. We have originally defined this cloud with 16 points 
surrounded by solid line. Within the dashed line, more 11 spectra are included to the left-top direction. 
Those 11 extra spectra, having irregular shapes in the same Vlsr component toward the peak position, 
may be part of this cloud. The dotted line cuts out only 4 spectra toward the central part. Then we have 
also calculated Leo and M^ii for both larger and smaller clouds. 
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Fig. 12. For eleven clouds, we have calculated Leo and Mvir at three different threshold values of CO 
integrated intensity in a same manner in Figure 1 1 , and all 33 points are plotted on the Lco--^vir diagram, 
compared with the Warp clouds (crosses), same as Figure 4. Thick line shows the least-squares fit for 55 
inner clouds identified by the NANTEN data. 
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